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ABSTRACT 

We investigate the Kozai mechanism in the context of near-Earth binaries and the Sun. The Kozai 
effect can lead to changes in eccentricity and inclination of the binary orbit, but it can be weakened 
or completely suppressed by other sources of pericenter precession, such as the oblateness of the 
primary body. Through numerical integrations including primary oblateness and 3 bodies (the two 
binary components and the Sun), we show that Kozai cycles cannot occur for the closely-separated 
near-Earth binaries in our sample. We demonstrate that this is due to pericenter precession around 
the oblate primary, even for very small oblateness values. Since the majority of observed near-Earth 
binaries are not well-separated, we predict that Kozai cycles do not play an important role in the orbital 
evolution of most near-Earth binaries. For a hypothetical wide binary modeled after 1998 ST27, the 
separation is large at 16 primary radii and so the orbital effects of primary oblateness are lessened. 
For this wide binary, we illustrate the possible excursions in eccentricity and inclination due to Kozai 
cycles as well as depict stable orientations for the binary's orbital plane. Unstable orientations lead 
to collisions between binary components, and we suggest that the Kozai effect acting in wide binaries 
may be a route to the formation of near-Earth contact binaries. 

Subject headings: minor planets, asteroids: general - minor planets, asteroids: individual (2002 CE26, 
2004 DC, 2003 YTl, Didymos, 1991 VH) 



1. INTRODUCTION 

About 15% of near-Earth asteroids (NEAs) larger than 
approximately 20 meters in diameter are in a binary 
configuration (Margot et al.l 120021 : iPravec et all I2006D . 
These NEA binaries have d ynamical lifet i mes o n the or- 
der of a few million years (|Bottke et al.l [20021 ) and are 
subject to a variety of perturbations, including close scat- 
tering encounters with terrestrial planets, radiative ef- 
fects (called BYORP) from the Sun, tidal torques, and 
the Kozai mechanism. Such perturbative effects may 
change the orbital energy and angular momentum of the 
system, which can influence the binary's orbital elements, 
including semi-major axis, eccentricity, and inclination. 
We briefly discuss each of these perturbations next; their 
relevance in explaining obser ved spin and orb i tal dat a of 
NEA binaries is discussed in iFang k MargotI ()2012bf ) . 

The effect and f requency of clos e plane tary encoun- 
ters are studied bv iFang fc Margo'd (|2012al ). who found 
that approaches (<10 Earth radii) with Earth can oc- 
cur for most observed NEA binaries on 1 — 10 million- 
year timescales. The radiative BYORP effect, while 
not observationally verified to date, is only relevant 
for satellites with spin-orbit synchronization. This ef- 
fect has been theoretically found to be capable of mod- 
ifying an NEA binary's semi-major axis and eccen- 
tricity on fast (~10^ years) timescales ()Cuk &: Burnd 



tion as well as niodify eccentricity (i.e. iGoldreich 119631: 
Gqldreich fc Saril 120091 ). The Kozai mechanism ( Kozail 
11962), in the context of NEA binaries, has not been fully 
studied and is the focus of this paper. 

The Kozai mechanism is very relevant for many astro- 
physical triple systems; examples include its influence 
on the stability for irregular Jovian sat ellites with 
high inclinations (e.g. INesvornv et al.l 120031). main bel t 
and trans- Neptunian binaries ( Perets fc Naod 2009D . 
asteroids and comets due to Jupiter fe.g. iKozai 119621 : 
iThomas & Morbidelli 1996), b inary stars with distant 
companions (|Harrington 1968; Mazch & Shaham 1971 
Kiseleva et al.l 19981: Eggleton fc Kiseleva-Eggleto 
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Fabrvckv fc Tremainel I2007D. 



outly i ng perturbers 



fc Kisseleva-Eggletonl 



. , Mazeh et al 
1997^ 



Eop 

extrasolar pl anets with 
^97; Inn anen et al.l 



1997; 'Holman et al."1997; 'Tremaine fc Zakamska"2004; 
Takeda fc Rasio 2005; Fabrvckv fc Tremaine ,2007 : 
20111: iLithwick fc Naod 120111: iNaoz et al 



Katz et al 

2011 ) , and binary supermassive black holes (jBlaes et al 



mm . 



20051 : ICuai2007l: IGoldreich fc Sarill2q09l : ICuk fc Nesvorn\^ 
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201C : iMcMahon fc ScheeresI l2010allbl : ISteinberg fc Sar 
201 It ). Tidal evolut ion of NEA binar i es ha s been pre- 
viously studied (see iTavlor fc MargotI ()2011l ) and refer- 
ences therein), and can cause spin-orbital synchroniza- 
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In this paper, we consider the Kozai mechanism in 
terms of the following triple system: an NEA binary, con- 
sisting of a massive primary and less ma ssive secondary, 
and th e Sun as an outer perturber (as in lPerets fc Naozl 
()2009() ). The Kozai mechanism is a secular effect, i.e. 
the effect occurs on timescales that exceed the orbital 
periods. Analytical computation of secular effects tra- 
ditionally involve averaging quantities over a complete 
orbital cycle. Under the Kozai effect, the secondary's 
orbit with respect to the primary will undergo coupled 
changes in eccentricity and inclination. We define the in- 
clination as the relative inclination between the binary's 
mutual orbit and the binary's heliocentric orbit. For a 
circular heliocentric orbit, the coupled oscillations in ec- 
centricity e and inclination i will conserve the quantity 



Table 1 

Sample of Near-Earth Binaries 



System 


Rp 


Mp 


Rs 


Ms 


n 


a 


e 


^Hill 




a© 


60 




(km) 


(kg) 


(km) 


(kg) 


(dog d-i) 


(km) 




(km) 


(deg d-i) 


(AU) 




(276049) 2002 CE26^ 


1.75 


2.17 X 10^^ 


0.150 


1.37 X lO^" 


554.34 


4.87 


0.025 


514 


0.30 


2.23 


0.56 


2004 DC^ 


0.17 


3.57 X lO^o 


0.030 


1.96 X 10* 


372.93 


0.75 


0.30 


44 


0.47 


1.63 


0.40 


(164121) 2003 YT1_^ 


0.55 


1.89 X 10^2 


0.105 


1.32 X 10^" 


226.73 


3.93 


0.18 


113 


0.84 


1.11 


0.29 


(65803) Didymos^ 


0.40 


5.24 X 10" 


0.075 


3.45 X 10^ 


724.38 


1.18 


0.04 


109 


0.47 


1.64 


0.38 


(35107) 1991 VH^ 


0.60 


1.40 X 10^2 


0.240 


8.93 X 10^° 


265.07 


3.26 


0.06 


105 


0.81 


1.14 


0.14 


Hypothetical Wide Binary 


0.42 


7.73 X 10" 


0.050 


1.51 X 10^ 


65.46 


6.66 


0.3 


62 


1.33 


0.82 


0.53 



This table consists of a subset of well-characterized NEA binaries whose uncertainty region for the orbital plane orientation includes 
the Kozai-acting regime (39.2° < i < 140.8°) for initially circular binaries. In addition, we list a hypothetical wide binary modeled after 
1998 ST27, whose actual physical and orbital properties are not well known. For all entries, we list adopted values for primary radius 
Rp, primary mass Mp, secondary radius Rs, secondary mass Ms, binary mean motion n, binary semi-major axis a, binary eccentricity 
e, the Hill radius i?Hill (beyond which the binary's components would be primarily orbiting the Sun instead of each other), heliocentric 
m ean motion np), heliocen tric semi-major axis oq, and heliocentric eccentricity sq. Binary parameters and uncertainties can be found 
in Fang & Margot ( 2012b'), and heliocentric parameters and uncertainties can be found in the JPL Small Body Database. 

= Shcpard et al. ( 2006,) 

t.Taylor ct al. ( 200^ 

Wolan ct al. (2004) 

'^Bcnner ct al. (201(J) 

^Margot et al. ( 200#) : rPravec et al.l ||2006| ') 



VT— f?cosi such that for satelhtes in prograde orbits 
{0° < i < 90°), peaks in eccentricity correspond to min- 
ima in incHnation, and vice versa. Satelhtes in retrograde 
orbits (90° < i < 180°) wih have eccentricity and incli- 
nation oscillate in the same direction. In the absence of 
other perturbations, an initially circular binary will un- 
dergo large Kozai cycles if a critical inclination i is met; 
39.2° < i < 140.8°. Initially eccentric binaries will un- 
dergo Kozai cycles over a wider inclination range. For an 
initially circular binary, the maximum eccentricity Cmax 
that can be induced by the Kozai mechanism is 



- C0S2 iinit 



(1) 



wher e ii^a is the initial inclination (i.e. Ilnnanen et alj 
119971 ). If Zinit ~ 90°, then the binary's eccentricity will 
grow to ~1 during a single Kozai oscillation. The ap- 
proximate os cillation timescale or Kozai period Pk is on 
the order of (|Kiseleva et al.lll998[ ) 



Pk = — ^(1 



Mr, +Ms+M(: 



(2) 



where the binary's mutual orbit has period P, the pri- 
mary's mass is Mp, the secondary's mass is Mg, the 
Sun's mass is Mq, the heliocentric orbital period is P©, 
and the heliocentric eccentricity is e©. The Kozai effect 
due to the Sun for binaries in the Solar System is much 
more pronounced for NBAs than for main belt and trans- 
Neptunian objects because Pk varies as the square of the 
heliocentric orbital period P©. 

Among the popu lation of well-cha r acteriz ed NEA sys- 
tems compiled by iFang fc MargotI (|2012b[ ). we list in 
Table [T] only binary systems for which the uncertainty 
region in orbital plane orientation includes the Kozai 
regime (39.2° < i < 140.8°) for initially circular bina- 
ries. This criterion rules out 2 triple systems, 2001 SN263 
and 1994 CC, and 2 binary systems, 2000 DP107 and 
1999 KW4 , for which the Kozai effect will not operate 
(|Fang fc Margot .2012b.) . We also include a hypotheti- 



cal wide binary modeled after 1998 ST27, which is the 
widest NEA binary known so far but unfortunately has 
binary orbital and physical parameters that are not well- 
determined. For the hypothetical binary's heliocentric 
orbital elements listed in TablelU we use 1998 ST27's ac- 
tual heliocentric orbital elements, which are well-known. 

The goal of this study is to investigate the relevance 
and effect of Kozai cycles in the presence of other mod- 
ulating perturbations. The outline of this paper is as 
follows. In Section [21 we describe other pcrturbations- 
the primary's oblateness, additional satellites, and tidal 
effects-that may damp Kozai cycles and we find that pri- 
mary oblateness is dominant. In Section [3l we perform 
numerical investigations of the Kozai effect modulated by 
primary oblateness for all binaries in our sample, and de- 
termine their excursions in eccentricity and inclination. 
In Sectional we discuss the implications of this work, in- 
cluding the role of Kozai cycles in the evolution of most 
NEA binaries, constraints on orbital plane orientations, 
how Kozai-induced instabilities end in collisions, and the 
formation of contact binaries. We briefly summarize this 
study in Section [31 

2. PERTURBATIONS THAT AFFECT KOZAI CYCLES 

Here, we consider effects that can weaken or completely 
suppress Kozai cycles. The Kozai effect causes oscil- 
lations in eccentricity and inclination, and systems in 
Kozai resonance can exhibit libration of the argument of 
pericenter. 

Given that this effect is caused by the interaction be- 
tween the shape of the binary's mutual orbit and weak 
solar tides, Kozai cycles can be easily suppressed by other 
weak perturbations that contribute to apsidal (pericen- 
ter) precession in the binary. If the argument of peri- 
center precesses too fast, libration of the argument of 
pericenter is no longer possible, which inhibits the Kozai 
process. 

In the following subsections, we consider contributions 
to pericenter precession from 3 effects: the primary's 
non-spherical shape leading to a non-uniform gravita- 
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tional field with a quadrupole moment (J2) representing 
the degree of oblateness, the presence of additional, unde- 
tected satellites, and tidal bulges due to both the primary 
and the secondary. These effects can potentially weaken 
or suppress Kozai cycles, and we seek to determine their 
relative strengths. 

2.1. Primary Oblateness (J2) 

We consider the contribution to pericenter precession 
caused by primary oblateness. Many NEAs have non- 
spherical shapes, and the leve l of oblateness can be de- 
scribed by a coefficient, J2, as ([Murray fc Dermottlll999D 



J2 - 



C-{A + B)/2 C-A 



MpRl 



(3) 



where A, B, and C are the primary's moments of inertia. 
Mp is the primary's mass and Rp is the primary's equa- 
torial radius. The approximation is valid when Ap^ B. 

The J2 coefficient is an indirectly observable quantity 
that can be detected by its non-Keplerian effects induced 
on orbiting satellites. Such oblateness-induced preces- 
sion on the orbits of satellites can be examine d by the 
rate of change in the argument of pericenter w (|Valladol 
[2OOI : 



duj 
dt J2 



2(1 



= 2^2 



Rr. 



^11- cos^ / - - 
a I \2 2 



(4) 



where n is the binary's mean motion, e is the eccentric- 
ity, Rp is the primary's radius, a is the semi-major axis, 
and / is the binary's orbital inclination relative to the 
primary's equator. The J2 effect is more relevant for as- 
teroid binaries than their trans-Neptunian counterparts 
because asteroid binaries tend to be separated by several 
primary radii and trans-Neptunian binaries are typically 
much wider. The orbits of close-in satellites will be more 
perturbed by an oblate primary than the orbits of distant 
satellites, since Equation 2] shows that precession due to 
J2 varies inversely as distance squared. 

The primary's J2 value for most NEA binaries is un- 
known; we calculate a range of values for apsidal preces- 
sion due to J2 coefficients ranging from 0.001 to 0.1. A 
few well-characterized NEA systems have known shapes 
and J2 values, including 1999 KW4 (lOstro et al.ll2006D 
with a primary J2 of ~0.06 and 1994 CC (jBrozovic et al.l 
[20Tl[) with a primary J2 of ~0.01. In our calculations, 
we assume that the satellite is orbiting in the equato- 
rial plane of the primary (support for this assumption 
comes from the generally accepted rotational fission for- 
mation mod el: iMargot et al.l 120021: iPravec et all 120061: 
IWalsh et al.ll2008D . and we use nominal values of their 
current separations and eccentricities. These numbers 
are presented in Table [2l 

2.2. Additional Satellite 

We examine if the presence of an additional satellite 
(presumably undetected) in the asteroid system can sup- 
press Kozai oscillations by causing the known satellite's 
orbit to precess. It is known that the presence of larger, 
detectable satellites can easily suppress K ozai cycles, 
as shown through numerical simulations by iFang et al.l 
(pOll.1 for NEA triple systems 2001 SN263 and 1994 CC 



and as discussed bv lRagozzine fc Brownl (|2009| ) for trans- 
Neptunian triple Haumea. Here, we investigate if a small, 
unobserved satellite can also damp Kozai cycles. 

Since most NEA binaries are discovered by planetary 
radar, we adopt a fairly common value of the spatial res- 
olution in radar images (15 m) as the typical radius of 
the additional satellite in our study. This is adopted for 
convenience and does not represent the finest spatial res- 
olution available nor the radar detectability threshold. 
If a small, undetectable satellite can damp Kozai cycles, 
then it may be responsible for the survival of NEA sys- 
tems that would otherwise undergo high oscillations in 
eccentricity and inclination. 

In the presence of an additional satellite, both satel- 
lites will undergo time-averaged, secular changes in their 
orbital elements. For a coplanar system , their argument 
of pericenter rates are (jMardlingj|2007fl 



dm si 

dt sat 



M. 



:nsi 



1-'- 
4 



s2 



Mp 

as2 



as2 
es2 
esi 



a-ei2r'' (5) 

C0s(Wsi - UJs2) 



1 



p2 
^s2 



dWs2 
dt sat 



r"-s2 



1 - 



5 / 

0-32 



asi 

0,32 

es2 



(1 

(1- 



4e 



(1 



(6) 

2 ^ COs(Wsi - LUs2) 



32_ 

32) 



and these equations are given to fourth power in 031/ as2 
and first order in esi. The subscripts si and s2 represent 
the inner and outer satellites, respectively, and the sub- 
script p is for the primary. The equations also include 
the mean motion n, mass M, semi-major axis a, and ec- 
centricity e. These equations are valid for M31 « Mp, 
but there are no restrictions on Ms2. 

We calculate the apsidal rates for each binary in Ta- 
ble [l] assuming a coplanar system. For each binary, the 
additional satellite is given a radius of 15 meters and a 
typical rubble pile density of 2 g cm~^ with a circular or- 
bit. For all binaries except the hypothetical wide binary 
modeled after 1998 ST27, we treat the observed satel- 
lite as the inner satellite and the additional, undetected 
satellite as the outer satellite with a semi-major axis of 
15 primary radii (which is a typical separation for the 
outer satellite in an NEA triple, based on a sample of 2 
known NEA triples). For the hypothetical wide binary, 
whose satellite is located at 6.66 km or 16 primary radii, 
we treat this satellite as the outer satellite and the un- 
detected satellite as the inner satellite with a semi-major 
axis of 4 primary radii. In all cases, we use the observed 
satellite's current separation and eccentricity. Our cal- 
culated apsidal precession rates are shown in Table [21 
which represent the pericenter precession due to a 30 m 
diameter satellite. 

2.3. Tidal Bulges 

For completeness we investigate tidal bulges raised on 
both the primary and secondary, although we anticipate 
their contribution to pericenter precession to be small . 
The argument of pericenter rate is (e.g. iSternd 119391 : 
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Table 2 

Pericenter Precession Rates 



System 


Pericenter 


Precession 


(deg d ^ 


)■■ 












CS 


2002 CE26 


0.22- 


-22 


<3.5e-6 


5.0C-7 - 


3.1e-5 


2004 DC 


0.069- 


-6.9 


<5.6e-3 


5.7e-8 - 


8.2e-5 


2003 YTl 


0.014- 


-1.4 


<2.8e-4 


8.6e-9 - 


l.Oe-6 


Didymos 


0.25- 


-25 


<2.2e-4 


1.3e-6 - 


3.0e-4 


1991 VH 


0.027- 


-2.7 


<1.9e-4 


3.3e-7 - 


8.3e-6 


HWB* 


0.00094-0. 


094 


<1.4e-4 


1.5e-ll- 


6.7e-9 



Rates for the argument of pericenter u) are calculated for 
each binary and for three main sources of pericenter pre- 
cession: primary's J2, an additional satellite, and tidal 
bulges raised on the primary and secondary. A range of 
rates is given for the effects of oblateness, correspond- 
ing to J2 values from 0.001 to 0.1. A range of rates is 
given for the effects of tides, corresponding to two dif- 
ferent tidal Love number m odels (Goldrcicli & Sari 20Q3; 
IJacobson &: Sch eeresI 120111). Precession duo to an addi- 
tional satellite assumes a size of 30 m. 

*HWB = hypothetical wide binary modeled after 
1998 ST27 

iFabrvckv fc Tremaind[2007t iBatvgin fc Laughlidl2011[ ) 



doj 

dt tides 



15 
7 

16 



8 + 12e^ 



(l-e2)5 



(7) 



Rr 



[Ms 



where the first term corresponds to the primary and the 
second term corresponds to the secondary. The equation 
also includes the mean motion n, eccentricity e, radius 
i?, semi-major axis a, mass M, and tidal Love number k. 
Subscripts p and s stand for the primary and secondary, 
respectively. 

The tidal Love number is a poorly known quantity for 
small rubble pile asteroids, and there are currently two 
different rubble pile models that descri be the Love num- 
ber's d ependence on size (or radius R) . iGoldreich fc Saril 
(|200t^ give the relat i on fcr ubbic ~ 10~^ km), and 

IJacobson fc ScheeresI (|201lD find that fcrubWo ~ 2.5x 10~^ 
(1 km/i?). Using both tidal Love number models and 
the binaries' current separations and eccentricities, we 
calculate a range of possible apsidal precession rates due 
to tidal bulges for each NEA binary in Table [21 

Examination of the values in Table [2] indicates that 
even the smallest amount of primary oblateness (J2 of 
0.001) will dominate over all other perturbations; the 
contribution from each perturbation is quantified in Ta- 
ble [2] for all NEA binaries in our sample. Therefore, 
the numerical integrations described in the next section 
only include the effect of primary oblateness out of the 3 
sources of pericenter precession considered here. 

3. NUMERICAL INVESTIGATION 

Through numerical simulations, we explore the excur- 
sion in a binary's eccentricity and inclination due to 3- 
body effects, such as Kozai cycles induced by the Sun, 
and we include perturbations due to primary oblate- 
ness. We use a Bulirsch-Stoer algorithm from Mercury 
(|Chamberslll999 l) and our system is composed of 3 bod- 
ies: the Sun and primary and secondary components of 



the binary. For each binary in our investigation (Table 
[1]), our integration time covers at least 10 Kozai oscilla- 
tion periods with a timestep that is 1/50*^ of the binary's 
mutual orbital period. Initial conditions for each binary's 
heliocentric semi-major axis and eccentricity as well as 
starting values for the masses, separations, and eccen- 
tricities of the binary components are taken from their 
known, observed values (Table |T]). 

For each binary in our sample, we perform an ensem- 
ble of simulations. In all simulations we assume that 
the binary's mutual orbit is in the primary's equatorial 
plane. We sample a range of J2 values to approximate 
the primary's non-spherical shape, a range of inclinations 
between the binary's mutual and heliocentric orbits, and 
a range of values for the binary's argument of pericenter; 
none of these parameters are known for the binaries in 
our sample. Our choice of J2 values includes 0.001, 0.005, 
0.01, 0.05, and 0.1. Curre ntly, the best-cha racterized 
NEA binary is 1999 KW4 (jOstro et all l2006t ) with de- 
tailed shape and orbital information, and its primary has 
a J2 value of ~0.06. Anot her NEA system with known 
primary shape is 1994 CC ( !Brozovic et al .''2011). with a 
primary J2 of ~0.01. Our range of inclinations spans an- 
gles in the strongest Kozai-operating regime of 40°— 90° 
with 10° increments; these inclinations represent satel- 
lites in prograde orbits. Identical behavior, but mirrored 
across 90°, would be seen for inclinations greater than 
90° representing retrograde orbits. We sample arguments 
of pericenter from 0° to 90° with 30° increments, as sim- 
ilar behavior is repeated for each 90° quadrant. 

Our findings are as follows. Due to a non-zero J2, 
binaries 2002 CE26, Didymos, 2004 DC, and 1991 VH 
do not exhibit any expected Kozai behavior such as li- 
bration of the argument of pericenter or large excur- 
sions in eccentricity and inclination. Instead, the argu- 
ment of pericenter circulates quickly and there are rapid, 
small-amplitude oscillations in eccentricity and inclina- 
tion. Typical behavior for these binaries that do not 
exhibit expected Kozai cycles are illustrated in Figure 
[TJs left panel. This figure displays the results from all 
numerical simulations for 2002 CE26 by showing the ec- 
centricity and inclination excursions. Moreover, if Kozai 
cycles are present, then in the absence of other perturbers 
we expect an initially circular binary to have an eccen- 
tricity increase (see Equation[T]) of ~^0.15 for jjnit of 40°, 
~0.56 for iinit of 50°, -0.76 for ijnit of 60°, -0.90 for ii„it 
of 70°, ~0.97 for iinit of 80°, and 1 for ijnit of 90°. Com- 
parison between these values and Figure [l]s left panel 
for 2002 CE26, an initially near-circular binary, provides 
additional evidence that the expected Kozai cycles are 
not present. For binary 2003 YTl, there are signs of 
Kozai cycles only with a J2 as low as 0.001, where there 
are significant excursions in eccentricity and inclination; 
larger J2 values suppress any Kozai oscillations. 

Kozai cycles may induce collisional disruptions be- 
tween the primary and secondary if the eccentricity grows 
large enough that the pericenter approaches one primary 
radius. This is more likely to occur at high inclinations 
(Equation [J). For binaries 2002 CE26, 2004 DC, Didy- 
mos, 1991 VH, and 2003 YTl, collisional disruptions are 
expected at high inclinations yet not observed in any of 
our simulations for the entire explored range of J2 val- 
ues. None of the simulations for each of these binaries 
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show any disruptions (i.e. collisions or ejections) during 
the nominal integration time, contrary to expected Kozai 
behavior that leads to high eccentricities within one cy- 
cle. The lack of disruptions in non-zero J2 cases provides 
further confirmation that expected Kozai behavior is not 
present when we include the effects of primary oblate- 
ness. 

For the hypothetical wide binary modeled after 
1998 ST27, Kozai cycles are present with the complete 
range of J2 values sampled here. In the case of high J2 
values, the system exhibits wide, overlapping excursions 
in eccentricity and inclination (Figure [TJs right panel). 
These Kozai cycles are modulated by a non-zero J2', 
Figure [2] illustrates a comparison of pure Kozai cycles 
vs. Kozai cycles under the influence of J2. This hypo- 
thetical binary exhibits Kozai behavior as opposed to 
the other binaries in our sample because it is the most 
widely-separated system with its satellite stationed at a 
relatively far separation of 16 primary radii. Therefore, 
solar perturbations are stronger than the effects induced 
by primary oblateness. For this hypothetical wide binary 
only, we continue integrating its ensemble of numerical 
simulations to 10^ years to obtain more accurate disrup- 
tion statistics. We find that Kozai-induced disruptions 
occur in ^AA% of all cases examined here, with colli- 
sions and not ejections as the only observed disruption 
outcomes. See Section |4] for discussion regarding such 
instabilities. 

These numerical integrations suggest that the evolu- 
tion of most observed NEA binaries are dominated by 
primary oblateness rather than the Kozai effect. This 
conclusion is also supported by an analytical compari- 
son bet ween the two effects. Fo r all well-characterized 
binaries, iFang fc Margotl (|2012bl ) calculated the critical 
semi-major axis separating the influenc e regions of pri- 
mary oblateness and solar dynamics as ([Nicholson et al.l 
[200^ : 

acHt=(2J2^i?X)'^' (8) 

where J2 represents the extent of primary oblateness, Mp 
is the primary's mass, Mq is the Sun's mass, Rp is the 
prima ry's radius, and is th e heliocentric semi-major 
axis. IFang fc Margotl (|2012bD found that the range of 
allowable Odit values corresponding to J2 values ranging 
from 0.001 to 0.1 were greater than observed semi-major 
axes for all well-characterized NEA binaries; they are 
oblateness-dominated. We also calculate Odit for the hy- 
pothetical wide binary modeled after 1998 ST27. For 
this hypothetical binary, Odit ranges from 3.02 — 7.60 
km and so could be larger or smaller than the adopted 
semi-major axis of 6.66 km. Therefore, this hypothetical 
binary could be in the oblateness-dominated or Kozai- 
dominated regime. This analysis shows general consis- 
tency with the results from numerical simulations; bina- 
ries predicted to be oblateness-dominated by Equation 
[8] do not show expected Kozai oscillations in numeri- 
cal simulations, and the hypothetical wide binary pre- 
dicted to be either oblateness or solar dominated does 
show strong signs of Kozai cycles modulated by J2 . Ac- 
cordingly, tightly-bound binaries (those that exhibit no 
signs of the Kozai effect) are strongly affected by primary 
oblateness and loosely-bound binaries (showing Kozai os- 



cillations) are less perturbed by the far-away primary's 
oblateness. 

4. DISCUSSION 

In this section, we discuss each of the following in turn: 
the relevance of Kozai cycles for observed NEA binaries, 
orbit orientation constraints for wide binaries that can 
undergo Kozai oscillations, collisions as the only disrup- 
tion outcome due to the Kozai effect for typical NEA 
binaries, and Kozai cycles as a possible route to the for- 
mation of contact binaries. 

4.1. General Population of Observed NEA Binaries 

We extrapolate from our simulation results to the 
observed NEA binary population. Numerical integra- 
tions in Section [3] demonstrated that 4 out of 5 well- 
characterized NEA binaries in our sample (Table [Ij 
showed no signs of Kozai cycles and the exception, 
2003 YTl, only showed evidence of Kozai cycles when 
the primary had a very minimal level of oblateness 
(J2 =0.001). All 5 of these well-characterized binaries 
have semi-major axes less than 8 primary radii. Only 
the hypothetical wide binary (at 16 primary radii) mod- 
eled after 1998 ST27 exhibited signs of Kozai cycles at 
a range of J2 values from 0.001—0.1. The effects of pri- 
mary oblateness are strongest for close-in binaries (peri- 
center precession due to J2 increases as the semi-major 
axis squared; Equation 14]) , where an oblate primary can 
cause the satellite's orbital plane to precess fast enough 
to thwart any Kozai oscillations. 

If we cons ider all well-characteriz ed NEA binaries as 
compiled bv lFang fc Margotl ()2012b[ ). none of them have 
separations greater than 8 primary radii. In fact, the 
semi-major axes of all of these well-characterized NEA 
binaries are within the range of semi-major axes of the 
binaries in our sample, which did not show Kozai cycles 
at a range of plausible J2 values. If we consider this list 
of well-characterized NEA binaries to be representative 
of the observed population of NEA binaries, then typical 
NEA binaries do not appear to be affected by Kozai cy- 
cles in their orbital evolution. We point out a mild selec- 
tion effect in that most well-characterized NEA binaries 
have made close enough approaches to Earth to be de- 
tected by radar, and this subset of the population may 
have fewer wide binaries than the general population. 
One could argue that we mostly observe binaries immune 
from Kozai effects because the Kozai-acting binaries have 
disrupted on short Kozai timescales (see following sub- 
sections on disruptions and contact binary formation). 
However, our simulations show that the requirements on 
primary oblateness and component separation are quite 
stringent, and we conclude that Kozai cycles are unlikely 
to be an important effect in most NEA binaries. 

4.2. Constraints on Orbit Pole Orientations 

Orbit pole orientations can be constrained for binaries 
that undergo Kozai cycles. In Section [3l we explored a 
hypothetical wide binary modeled after 1998 ST27 that 
exhibited signs of Kozai oscillations and disruptions, and 
in this subsection we illustrate this binary's stable orbit 
pole orientations. Stable orbit poles are defined as ini- 
tial orientations that do not end in a collision after 10^ 
years of numerical integrations; these results are based 
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Figure 1. Results from all numerical simulations for 2002 CE26 (left) and for a hypothetical wide binary modeled 
after 1998 ST27 (right) are shown in these stacked plots with panels representing different J2 values for the primary. 
The y-axis shows excursions in eccentricity e space and the x-axis illustrates the ranges of inclination i. In each panel, 
starting conditions include the observed separation and eccentricity as well as various values for the inclination and 
the argument of pericenter. 

on the ensemble of simulations performed in Section [3] 
as well as additional simulations to sample the full range 
of inclinations. The range of inclinations correspond- 
ing to stable binaries allows us to calculate the corre- 
sponding J2000 ecliptic coordinates (latitude (3, longi- 
tude A) of binary orbits for stable binaries. Since the 
hypothetical wide binary under consideration is modeled 
after 1998 ST27, we use 1998 ST27's well-known helio- 
centric orbit poles (ecliptic longitude of the ascending 
node=197.5842°, ecliptic inclination=21.05458°) for our 
calculations. 

In Figure [3l we map out constraints on the orbit orien- 
tations for this hypothetical wide binary, showing results 
for 3 different J2 values for the primary: 0, 0.005, and 
0.05. Color-coded lines (for different J2 values) separate 
the "Kozai stable" and "Kozai unstable" regions. Re- 
gions where at least half of simulations resulted in sta- 
ble binaries are called "Kozai stable," and regions where 
over half of simulations resulted in unstable binaries are 
called "Kozai unstable." This orbit orientation map pro- 
vides constraints on allowable orientations because bina- 
ries that will disrupt under the effect of Kozai are unlikely 
to be observed with orbit ecliptic coordinates in the un- 
stable regions of Figure |3l Most Kozai-induced disrup- 
tions occurred within 100 years, shorter than evolution- 
ary timescales due to tides, BYORP, and close planetary 
encounters. Accordingly, for this hypothetical wide bi- 
nary we can predict that its orbit does not lie in the 
region bounded by the color-coded J2 lines. This anal- 
ysis is performed here for this hypothetical wide binary 
and can similarly be applied to other binaries that can 
be affected by Kozai perturbations. We note that a dif- 
ferent binary (presumably with a different heliocentric 
orbital pole) would result in a different layout of stable 
and unstable zones (i.e. stability islands). 

Figure [3] also shows how different J2 values can affect 
the stability map, discussed here for the hypothetical 
wide binary. For discrete J2 values of 0, 0.001, 0.005, 
0.01, and 0.05, we find that an increasingly larger J2 
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Figure 2. Here is a comparison between pure Kozai cycles 
(J2 = 0) and modulated Kozai cycles (J2 = 0.01) for a 
hypothetical wide binary modeled after 1998 ST27. The 
starting conditions for this particular simulation are as 
follows: eccentricity is 0.3, inclination is 40°, and the 
argument of pericenter is 90°. 
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value results in an increase in the number of unstable 
orbit orientations in our simulations. These higher J2 
values cause collisions that can occur at a larger range of 
initial inclinations by increasing the maximum eccentric- 
ity acquired during a sequence of Kozai cycles. A larger 
maximum eccentricity therefore decreases the pericen- 
ter distance to allow for more frequent collisions (e.g. 
iNesvornv et al.|[2003[ ). We do not continue to observe 
this trend at high J2 values greater than ~0.05, where 
substantial primary oblateness prevents additional un- 
stable orientations from forming. 

Hypothetical Wide Binary: Orbit Pole Constraints 



^ 100 



D -50 



J, values: —0 —0.005 — 0.05 



Kozai stable 




100 200 300 

Initial Ecliptic Longitude X (deg) 

Figure 3. This stability map shows constraints on ini- 
tial orbit pole orientations for a hypothetical wide binary 
modeled after 1998 ST27. There are 2 regions: "Kozai 
stable" refers to regions where at least 50% of simula- 
tions resulted in stable binaries, and "Kozai unstable" 
refers to regions where more than 50% of simulations re- 
sulted in unstable binaries. Intermediate lines of varying 
colors (representing different J2 values) separate "Kozai 
stable" and "Kozai unstable" regions. Stable orbit ori- 
entations cover roughly ~75%, ~50%, and ^45% of the 
celestial sphere for J2 values of 0, 0.005, and 0.05, respec- 
tively. Stability is defined as binaries with no disruptions 
after 10^ years in the ensemble of numerical integrations 
(Section ig. 



4.3. Kozai-Induced Disruptions End in Collisions 

Here we discuss how disruptions due to the Kozai ef- 
fect will result in collisions for typical NEA binaries. We 
only consider potential cases where the Kozai effect oper- 
ates and can produce instabilities, and we assume these 
instabilities occur on much shorter timescales (on the 
order of tens to hundreds of years) compared to other 
perturbations such as planetary encounters, tidal evolu- 
tion, and BYORP. We define disruptions as dynamical 
instabilities due to the Kozai effect, whose outcomes in- 
clude (1) collisions between binary components and (2) 
ejections where the binary becomes unbound. The Kozai 
effect is capable of increasing a binary's eccentricity, and 
very high eccentricities can cause the binary to become 
unbound or cause a collision between the binary compo- 
nents. An ejection occurs when an orbit's apocenter Q, 
which is related to the semi-major axis a and the eccen- 
tricity e as Q = a(l -|- e), grows to a distance greater 



than the binary's Hill radius i?Hin- A collision occurs 
when an orbit's pericenter q, where q — a{l — e), is less 
than the primary radius Rp. Rmii and Rp are known 
quantities given in Table [T] Therefore, for given i?Hiii! 
Rp, and a, whether a gradually increasing eccentricity 
(i.e. due to the Kozai effect) causes a binary to undergo 
a collision or ejection is solely determined by the value 
of the eccentricity. 

For all NEA binaries studied here (Table [1]), values 
for their maximum Q (by assuming an e of 1) are much 
smaller than their Hill radii. In order for Q to be as 
large as the Hill radius, if we assume a maximum e of 1 
then any binary's semi-major axis needs to be at least 
half as large as the Hill radius in order for an ejection to 
occur. Thus, for any typical binary with a semi-major 
axis smaller than 0.5i?Hin, Kozai cycles will not cause 
ejections. On the other hand, collisions will occur at an 
eccentricity less than 1. Collisions will occur (by set- 
ting q = Rp) at the following e values for the binaries 
in our sample: 0.64 for 2002 CE26, 0.77 for 2004 DC, 
0.94 for the hypothetical wide binary, 0.86 for 2003 YTl, 
0.66 for Didymos, and 0.82 for 1991 VH. The specific 
case of the hypothetical wide binary is most interesting 
because we have already shown that the other binaries 
do not generally undergo Kozai cycles, let alone Kozai- 
induced instabilities. But for any wide NEA binary such 
as the hypothetical case examined here, if Kozai cycles 
are acting with sufficiently high inclinations to induce 
high eccentricities, collisions between binary components 
are the only possible instability outcome. The binary's 
pericenter will shrink below the primary radius before 
its apocenter can increase beyond the Hill radius and be- 
come unbound. Since collisions are the only instability 
outcome, this means that the Kozai effect cannot directly 
form asteroid pairs. It is possible that the Kozai effect 
may indirectly form asteroid pairs if a collision occurs but 
the components separate again and become unbound, or 
if the binary fortuitously makes a close planetary en- 
counter during the short Kozai timescale. 

For Kozai cycles in the absence of other perturbations, 
these collisions occur when the starting inclination re- 
sults in a maximum eccentricity (Equation [ij that is 
greater than the limiting eccentricities listed in the pre- 
vious paragraph. Thus, for these starting inclinations, 
collisions will occur over the course of one Kozai cycle 
and repeated cycles will not occur. These disruption 
timescales are fast; for the hypothetical binary modeled 
after 1998 ST27, our longer-term (10^ years) simulations 
show that from a starting eccentricity of 0.3, collisions 
occur ranging from ~3 years to >10^ years later with 
the majority of collisions occurring within 100 years. 

4.4. Formation of Contact Binaries 

We discuss contact binary formation from large- 
amplitude Kozai oscillations. For wide binaries such 
as 1998 ST27 that can potentially undergo Kozai cycles 
even in the presence of primary oblateness, high starting 
inclinations can lead to high eccentricities (Equation [IJ . 
During a single Kozai oscillation, these high eccentricities 
can drive the orbital pericenter distance to very low val- 
ues and can cause the binary components to collide (see 
previous section for discussion on collisions). If a collision 
does not occur, tidal friction may play an important role 
and decrease the semi-major axis as well as the eccentric- 
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ity. When the pericenter distance decreases due to high 
eccentricities, tides can more efficiently circularize the or- 
bit by dissip ating orbit al energy dur ing each pericenter 
passa ge (e.g. lPerets fc N aoz 2009; Fa brvckv fc Tremaind 
l2007f l. This leads to a circular orbit and a smaller semi- 
major axis. 

If Kozai cycles and possibly tidal effects can cause bi- 
nary components to come into contact with each other, 
this may be a mechanism to create near-Earth contact 
binaries, which consti tute ^10% of N EAs larger than 
200 m in diameter ( Be nner et al J 120061 ). If this process 
occurs, we would expect the observed contact binaries 
to have high obliquities preferentially between ~40° and 
~140°, assuming no significant obliquity evolution has 
occurred after formation of the contact binary. Obliq- 
uity is defined here as the inclination between the con- 
tact binary's rotational angular momentum vector and 
the vector normal to its heliocentric orbital plane. This 
obliquity angle is equivalent to the inclination angle de- 
fined in Section[l]and used throughout the paper; we as- 
sume that the direction of an initially-detached binary's 
orbital angular momentum vector is the same as the di- 
rection of a then-collided contact binary's rotational an- 
gular momentum vector. 

Observed nea r-Earth contact bin aries include Castali a 
(iHudson fc Ost ro 1994) Bac chus (iBenner et al.l I1999D 
Mith ra (|Brozovic et al.ll2010D. 1996 HWl (|Magri et al.l 
120111 ). and Itokawa fe.g. lOstro et al.ll2004D : ah of which 
have high obliquities (Busch et al., in prep.). In the 
inner main belt, small binaries are observed with obliq- 
uities concentrated toward ° and 180° and the re is a 
lack of obhquities near 90° (jPravec et all 120 111) . The 
concentration of binaries with low inclinations and the 
concentration of contact binaries near high inclinations 
hints that the Kozai effect ma y be effective a t disru pting 
high-inclination binaries, but iPravec et al.l (|2011l ) have 
shown that the Kozai effect cannot completely explain 
the observed concentration of binaries with obliquities 
near 0° and 180°. 

Although the Kozai effect is not responsible for the 
observed binary pole concentrations of small main belt 
binaries near 0° and 180°, it can be effective for more 
rarely observed NEA binaries with wide separations. Our 
simulations show that widely-separated binaries can es- 
cape the dominating influence of primary oblateness to 
be significantly affected by solar perturbations. As a 
result, we suggest that wide binaries with eccentric or 
highly-inclined orbits may become unstable under the 
effects of Kozai, leading to the formation of near-Earth 
contact binaries. Alternate theories of contact binary 
formation include low-velocity collisions between com- 
ponents in an unstable binary due to planetary encoun- 
ters or radi ative effects such as YORP and BYORP (e.g. 
iTavlor fc Margot.2011) . 

5. CONCLUSION 

We explored the effect of Kozai cycles caused by the 
Sun on a sample of NEA binaries. Kozai oscillations can 
be suppressed by significant sources of pericenter preces- 
sion; we identified 3 processes-primary oblateness (J2), 
presence of an additional satellite, and tides-that con- 
tribute to orbital precession. We determined that pri- 
mary oblateness is a dominant source of pericenter pre- 
cession for small near-Earth satellites. Accordingly, we 



performed numerical simulations to evaluate the strength 
of Kozai cycles for NEA binaries in our sample with the 
inclusion of the primary's J2- 

Our study showed that the binaries in our sample (with 
the exception of a hypothetical wide binary modeled af- 
ter 1998 ST27) do not undergo Kozai cycles due to their 
small component separations. This sa mple includes all 
well-ch aracterized binaries as defined bv lFang fc MargotI 
(|2012bl ) but does not include all observed binaries. Even 
the presence of a minimal J2 (0.001) prevented NEA bi- 
naries in our sample from exhibiting signs of Kozai cy- 
cles. Consequently, we conclude that the Kozai effect is 
not relevant in explaining the observed characteristics of 
typical, observed NEA binaries. However, we note that 
the Kozai effect may have shaped the observed popula- 
tion of binaries by eliminating those binaries with very 
low- J2 primaries or widely-separated components. 

For rarer observed cases of wide binaries, we studied a 
hypothetical wide binary modeled after 1998 ST27 whose 
large component separation indicated weaker J2 effects 
and visible Kozai cycles at a wide range of possible J2 
values. For this hypothetical wide binary, we were able 
to map out orbit pole orientations that are stable un- 
der Kozai effects. For cases of such wide binaries, the 
Kozai effect can lead to collisions between components. 
Accordingly, we suggest that the Kozai effect acting on 
widely-separated binaries may be a route to the forma- 
tion of near-Earth contact binaries. 

We thank Simon Porter and the referee, Matija Cuk, 
for useful discussions. This work was partially supported 
by NASA Planetary Astronomy grant NNX09AQ68G. 
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